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Abstract. We prove decidability of univariate real algebra extended 
with predicates for rational and integer powers, i.e., £ Q” and 

“r" € Z.” Our decision procedure combines computation over real alge¬ 
braic cells with the rational root theorem and witness construction via 
algebraic number density arguments. 


1 Introduction 

From the perspective of decidability, the reals stand in stark contrast to the 
rationals and integers. While the elementary arithmetical theories of the integers 
and rationals are undecidable, the corresponding theory of the reals is decidable 
and admits quantifier elimination. The immense utility real algebraic reasoning 
hnds within the mathematical sciences continues to motivate significant progress 
towards practical automatic proof procedures for the reals. 

However, in mathematical practice, we are often faced with problems involv¬ 
ing a combination of nonlinear statements over the reals, rationals and integers. 
Consider the existence and irrationality of \/2, expressed in a language with 
variables implicitly ranging over K: 

Bx(x > 0 A x^ = 2) A -'Bx(x sQAa:>0Aa;^=2) 

Though easy to prove by hand this sentence has never to our knowledge been 
placed within a broader decidable theory so that, e.g., the existence and irra¬ 
tionality of solutions to any univariate real algebra problem can be decided auto¬ 
matically. This \/2 example is relevant to the theorem proving community as its 
formalisation has been used as a benchmark for comparing proof assistants mi¬ 
lt would be useful if such proofs were fully automatic. 

In this paper, we prove decidability of univariate real algebra extended with 
predicates for rational and integer powers. This guarantees we can always decide 
sentences like the above, and many more besides. For example, the following 
conjectures are decided by our method in a fraction of a second: 

yx{x^ GZAx^^Z => X ^Q) 

3x{x'^ GQAa;^QAa;^-|-l> 20) 

Va;{a;^ ^ Q => x ^ Q) 

3x{x ^ Q A £ Z A Sx"^ -I- 2x -f 1 > 5 A 4x^ -I-1 < 2) 


2 Preliminaries 


We assume a basic grounding in commutative algebra. We do not however as¬ 
sume exposure to real algebraic geometry and give a high-level treatment of the 
relevant foundations. 

The theory of real closed fields (RCF) is Th{{R, , x, <, 0,1)), the collec¬ 
tion of all true sentences of the reals in the elementary language of ordered rings. 
RCF is complete, decidable and admits effective elimination of quantifiers [3]. 

A real algebraic number is a real number that is a root of a (non-zero) 
univariate polynomial with integer coefficients. The real algebraic numbers, 

^aig = {x S K I ^ 0 S Z[a;] s.t. p{x) = 0}, 

form a computable subfield (a computable sub-RCF) of R. Indeed, Raig embeds 
isomorphically into every RCF. The field operations of RaZg are performed on 
computable representations of field elements. The minimal polynomial of a G 
Raig is the unique monic p G Q[a:] of least degree s.t. p{a) = 0. The degree of an 
algebraic number is the degree of its minimal polynomial. 

An element a € Maig can be represented by two pieces of data: (i) a polyno¬ 
mial p(a::) G Zlx] s.t. p(a) = 0, and (ii) an identifier specifying which root of p(x) 
is denoted by a. A root-triple representation is often used where a is “pinned 
down” among the roots of p(x) by an interval with rational endpoints: 

(p(x) G Z[x],qi,q 2 G Q) s.t. p{a) = 0 A #{r G [< 71 , 92 ] | p{r) = 0} = 1. 

The process of root isolation is a key component of computing over RaZg- 
Given a polynomial p G Z[a:] with k unique real roots, root isolation computes 
a sequence of disjoint real intervals with rational endpoints Ji,... ,7^ s.t. each 
Ij contains precisely one real root of p. Much work has been done on efficient 
root isolation. Common approaches include those based on Sturm’s Theorem 
and Descartes’ Rule of Signs mMM- Sturm’s Theorem also plays a key role in 
computing the sign of a polynomial evaluated at a real algebraic number. 

Given representations of a,/3 G RaZg, there are two main approaches to per¬ 
forming the field operations, i.e., for computing representations of a~^, oi-\-/3, ajS, 
etc. Both approaches rely on root isolation. The first approach uses bivariate re¬ 
sultants to compute representation polynomials [T3] . The second approach uses a 
recursive representation of real algebraic numbers through an explicit treatment 
of field towers and does not require computing resultants PITT] . Computing o:” 
(which plays a key role in our decision procedure) can in general be done by 
repeated squaring, requiring on the order of log n real algebraic number multi¬ 
plications. More sophisticated methods for a" are also available [5]. 

The Intermediate Value Theorem (IVT) holds over every RCF. Armed with 
machinery for computing the sign of a polynomial p{x) G Z[x\ at a real algebraic 
point a G RaZg, the combination of IVT and root isolation can be used as the 
basis of a decision method for univariate real algebra. 






Consider 


ip{x) = 


k^ kn 


A V Qij 0) 

i=l j=l 


s.t. pij G Z[a;], Qij G {<,<,=,>, >}■ 


We can decide the satisfiability of p over R, i.e., whether or not 
(R,+,-, x,<,0,l) 1= 3x{ip{x)) 


in the following manner: 

— Let P = YlijPij G ^[x], the product of all polynomials appearing in ip. 

— Let ai < ... < ak & RaZg be all distinct real roots of P. 

— Then, the roots partition R into finitely many connected components: 

R = ]-oo,ai[U [oi] U]ai,Q! 2 [U ... U]afc_i,afc[U [au] U]afc,+oo[. 


— By IVT, the sign of each polynomial pij appearing in p is invariant over any 
component of the partitioning. 

— Thus, we can simply select one sample point from each component of the 
partitioning and obtain a sequence of 2fc + 1 real algebraic points S = 
{ri, ... ,rfc+i} C M^alg S.t. 

2/c+l 

(R,+,-, x,<,0,l) 1= 3x{p{x)) \f pin). 

i=l 

Now 3a:((/3(a:)) can be decided simply by evaluating p{x) at finitely many 
real algebraic points. The partitioning of R constructed above is called an 
algebraic decomposition induced by P (equivalently, by the polynomials Pij). 


3 Decision Procedure 

Our decision procedure extends the IVT-based method for univariate real algebra 
with means to handle predicates expressing the rationality and integrality of 
powers of the variable of the formula, i.e., (a:" G Q) and (x" G Z). As will be 
made clear (cf. Sec. [5]), the restriction of these predicates to powers of the variable 
is important: The method would fail if we allowed more general polynomials 
p{x) G Z[x] to appear in constraints of the form (p(x) G Q). 

Formally, we work over the univariate language of ordered rings C extended 
with infinitely many predicate symbols of one real variable: 

(x G Q), (x^ G Q), (x^ G Q), ... and (x G Z), (x^ G Z), (x^ G Z),... . 

We use Cqi to mean the resulting extended language and £q (resp. £z) to mean 
C extended only with the rationality (resp. integrality) predicates. 




We present a method to decide the satisfiability of quantifier-free £qz for¬ 
mulas over R. It suffices to consider formulas of the form 

(p{x) A r{x) 

where G £ is a formula of univariate real algebra and 

r = Fq A Fz 


s.t. 


and 


Fn = 


A(= 




A 


^2 


.W2(i) 


■ fca ki 

G Z) A /\(a 


^Z) 


Informed by the IVT-based method for univariate real algebra, we can re¬ 
duce this £qz decision problem to an even more restricted one. Crucial to this 
reduction is treating the connected components of an algebraic decomposition 
as “first class” objects, rather than only computing with single sample points 
selected from them. We call such components r-cells. 


Definition 1 (r-cell). An r-cell is a connected component o/R of one of the 
following four forms (with a,(3 G Raigj- (i) [o;], (ii) ]—oo,a[ s.t. a < 0, (Hi) 
]a, /3[ s.t. 0 < a < (3 or a < (3 <0, (iv)]a, -|-oo[ s.t. a > 0. 

Observe that the only r-cell containing zero is the singleton (type (i)) r-cell 
[0]. Note that r-cells of type (i) are 0-dimensional subsets of M while r-cells 
of types (ii)-(iv) are I-dimensional. We call these 0-cells and I-cells, resp. An 
algebraic decomposition can always be transformed into an r-eell decomposition 
by splitting any I-cell containing zero into three parts. 

Given <l>{x) = ip{x) A F{x), we must decide whether or not R contains any 
point X s.t. <P{x) holds. To do so, we will first compute an r-cell decomposition 
of R induced by the polynomials of (p. Let ci,..., be these r-cells. Then by 
IVT, the truth of p is invariant within each c^. Note, however, that the truth of 
F may vary over each c^. Let C be the result of filtering out all r-cells Ci that 
falsify p: 

C = {ci I 3r G Ci{p{r)), 1 <i < k}. 

This can be done by evaluating at a single sample point drawn from each Cj. 
If C = 0, then <P is clearly unsatisfiable over R. Otherwise, C is a non-empty 
collection of r-cells over which p is satisfied. To decide F, we need only to decide 
whether or not F is satisfied over any c G C. 

We present a method to do so. We first develop a method to decide rationality 
constraints over an r-cell. We then lift the method to handle general combinations 
of rationality and integrality constraints. 



3.1 Deciding rationality constraints 

Given a system of rationality constraints Fq and an r-cell c, we need a method to 
decide whether or not Iq is satisfied over c. To accomplish this, we will extract 
a system of degree constraints from Fq and give a method to decide if c contains 
a real algebraic number satisfying them. 

We must however take care of the following issue: If we prove there exists no 
algebraic real in c satisfying Tq, how do we know there exists no transcendental 
real in c satisfying Fq as well? That is, in the presence of rationality constraints, 
can we still transfer results from Raig to M as a whole? We answer this question 
in the affirmative by proving a suitable transfer principle (cf. Theorem [5]). 

It turns out we need essentially two methods for deciding Fq over c: One 
method for 0-cells and another for 1-cells. We begin with the 1-cell case. 

1-cells To construct our system of degree constraints, we shall utilise a funda¬ 
mental property relating the degree of a “binomial root” real algebraic number 
to the rationality of its powers. We employ a result on the density of real al¬ 
gebraic numbers to show that any consistent system of degree constraints gives 
rise to a real algebraic solution in a 1-cell. We then prove completeness of the 
method and a transfer principle enabling us to lift results from ]&aig to R. 

Lemma 1 (Minimal binomials). Let a G Ro/g s.t. a" S Q for some n G N. 
Then, the minimal polynomial for a over Q[x] is a binomial of the form x‘^ — q. 

Proof. Let k G N be the least power s.t. G Q. We shall prove that p{x) = 
x^ — G Q[a:] is the minimal polynomial for a. Assume p{x) is reducible over 
Q[a:]. Observe that p{x) = where f is a kth root of unity. As p{x) 

is reducible, it must have a nontrivial factor f{x) — ~ £ QN with 

m < k and Si G N. But then (a™ YYiLi C*‘) G Q» ^.nd since a is real, we must 
have a™ G Q. But m < k. Contradiction. Thus, as p{x) = x^ — is irreducible 
and monic, it is the minimal polynomial for a over Q[a:]. □ 

Lemma 2 (Binomial algebraic degree and divisibility). Let a G KaZg s.t. 
a is a root of some x^ — q G Q[a:]. Let n G N. Then, 

(a" G Q) deg{a) \ n. 

Proof. Let d = deg{a). (<G=) By Lemma [IJ a‘^ G Q. But, as d \ n, we have a" = 
(a^)* for some k G N. Thus, a" G Q. (=>) We use the method of infinite descent. 
Consider a" = g G Q. Then, x^ — q has a as a root, and thus d < n. Assume 
d\ n. It follows that d < n, gcd{d,n) = 1, q = and gcd{d,n — d) = 1. 

As G Q, we have d^~‘^ = ^ G Q- Note n — d < n. But then d^~‘^ G Q s.t. 
d \ n — d, and we can continue this process ad infinitum. Contradiction. □ 

Let c C R be a 1-cell and Fq a system of rationality constraints s.t. 

' ki /C2 

/\(a:“^« G Q) A /\(:r“'=« ^ Q) . 

. 2=1 2=1 


n 




To Fq, we associate a system of degree constraints 'D^Fiq) as follows: 


P(/q) 


ki k2 

/\{d I A /\{d\w2ii)) 


.i=i 


i=l 


Note that each Wj(i) is a concrete natural number. Thus, T>(/q) is a system of 
arithmetical constraints with a single free variable d. We shall prove that /q is 
satisfied over c iff T’(Tq) is consistent over N, i.e., iff 


3dG Ns.t. V{rq){d). 


We proceed in two steps. First, we prove that Fq is satisfied by a real algebraie 
number in c iff 'D{Fq) is satisfied over N. Next, we show that this result can be 
lifted to R as a whole, i.e., that Fq is satisfied over c (by any real, be it algebraic 
or transcendental) iff 'F>{Fq) is satisfied over N. 

These results elucidate a deep homogeneity of K. Intuitively, M is so saturated 
with real algebraic numbers that, given any open interval / C M, the only way I 
can fail to contain an algebraic number satisfying Fq is if the purely arithmetical 
facts induced by Fq (via Lemma [5]) are mutually inconsistent over N. Moreover, 
from the perspective of rationality constraints, transcendental elements cannot 
be distinguished from algebraic ones. To prove these results, we shall need to 
understand a bit about the density of real algebraic numbers of arbitrary degree. 

Lemma 3 (Density of ratios of primes). Given a < & G R, there exists 
^ G ]a, 6[ s.t. \p\ ^ |g| are both prime. 

Proof. A straightforward application of the Prime Number Theorem. 

Lemma 4 (Density of real algebraic numbers of degree n). Let a < 6 G R 

and n G N. Then, 3a G MaZg s.t. a < a < b and deg{a) = n and a^ G Q. 

Proof. We construct an irreducible p{x) = — q € Q[a;] s.t. a < < b. 

Then, a = will suffice. WLOG, assume a > 0. Let Q be a rational in 
]a, 6[. Let f : R’*' —>■ R 6e the nth-root function, i.e., f{r) = l/r. Consider 
Q" G Q. By continuity of f, 3e > 0 s.t. /(]Q" — e, Q" + e[) C ]a, 6[. For each 
rational q G JQ" — e, Q" + e[, we thus have a < f{q) < b with f{q) algebraic, as 
— g = 0. To prove the theorem, we must choose q s.t. deg{f{q)) = n. It 
suffices to find q G ]Q" — e, Q" + e[ s.t. p(x) = x'^ — q is irreducible over Q[a;]. By 
Lemma\^ we can choose g = — G JQ" — e,Q" + e[ s.t. gi ^ g 2 are both prime. 
By Eisenstein’s criterion, g 2 a:" — gi is irreducible over Q[x]. Thus, 2 :" — ^ is 

irreducible and a = completes the proof. □ 

With LemmaSlin hand, it is not hard to see that Fq is satisfied by a real algebraic 
number in a 1-cell c iff 'D{Fq) is satisfied over N. 

Theorem 1 (1-cell arithmetical reduction: algebraic case). Let Fq be a 

system of rationality constraints and c C R a 1-cell. Then, Fq is satisfiable over 
c by a real algebraic number iff'D(Fq) is satisfiable over N. 




Proof. (=>) Let a G (cflKaZg) satisfy Pq. Then, by Lemma\^ d = deg{a) satisfies 
X>(/q). (<^=) Let d € N satisfy 2?(/q). Then, by Lemma\^ any algebraic a € c 
s.t. deg(a) = d will satisfy Pq. But, by Lemma\^ such an a must exist in c. □ 

Thus, we have reduced the satisfiability of Pq by real algebraic numbers 
present in a 1-cell c to the satisfiability of 'D{rQ) over N. However, we must still 
attend to the possibility that Pq could be satisfied by a transcendental element 
in c without being satisfied by an algebraic element in c. Let us now prove 
that this scenario is impossible. In fact, we will prove this for both the 0 and 
1-dimensional cases. 


Theorem 2 (Rationality constraints transfer principle). Let Pq be a sys¬ 
tem of rationality constraints and c an r-cell. Then, it is impossible for Pq to he 
satisfied by a transcendental real in c without also being satisfied by an algebraic 
real in c. 


Proof. Let Pq = 


If c is a 0-cell, then 


c contains no transcendental elements, so the theorem holds. Consider c a 1- 
cell. We examine the structure of Pq. If ki > 0, i.e., Pq contains at least 
one positive rationality constraint, then Pq cannot he satisfied by any transcen¬ 
dental element, and the theorem holds. Thus, we are left to consider Pq = 
^ Q) Aj *5 satisfied by a transcendental element in c. Let 
m = max(w 2 {C),... ,W 2 (k 2 )). Then, Pq will be satisfied by any a G Maig s.t. 
deg(a) > m. But by Lemma^ c must contain an algebraic a s.t. deg(a) = m-|-l. 

□ 


In addition to giving us a complete method for deciding the satisfiability of 
systems of rationality constraints over 1-cells, the combination of Theorem 
and the completeness of the theory of real closed fields tells us something of a 
fundamental model-theoretic nature: 

Corollary 1 (Transfer principle for Cq). Given (p G Cq, 

(R, -I-, X , <, {x^ G Q)r!,eN, 0, 1) ^ {Ralg,+, X , <, {x^ G Q)neN, 0, 1) \= (j). 

That is, extending the language £ to include rationality constraints (£q) still 
guarantees a sound transfer of results from Raig to M. 

Finally, let us put the pieces together and prove our main theorem for 1-cells. 

Theorem 3 (1-cell arithmetical reduction: general case). Let Pq be a 

system of rationality constraints and c C R. a 1-cell. Then, Pq is satisfiable over 
c iff 'D{Pq) is satisfiable over N. 

Proof. Immediate by Theorem [7] and Theorem □ 

Thus, to decide if Pq is satisfied over a 1-cell c, we need only check the 
consistency of £’(£q) over N. It is easy to derive an algorithm for doing so. 
Consider T>{Pq) s.t. 

~ ki k2 

/\{d I wi{i)) A /\{d\w 2 ii)) . 

. 2=1 2=1 


V{Pq) 






If ki = 0, then d = max{w 2 {l), ■ ■ ■ ,^ 2 (^ 2 )) + 1 satisfies 'D{rq). If ^2 = 0, then 
d = 1 satisfies 2?(rQ). Finally, if fci >0 and ^2 > 0, then m = mm(wi(l),..., wi{ki 
gives us an upper bound on all d satisfying ^(rq). Thus, we need only search 
for such a d from 1 to m. For efficiency, we can augment this bounded search by 
various cheap sufficient conditions for recognising inconsistencies in 2 ?(/q). 


0-cells When deciding rationality constraints over r-cells of the form [a], we 
will need to decide, when given some j G N, whether or not a-' G Q. Recall that 
a root-triple for can be computed from a root-triple for a (cf. Sec. [5]). A key 
component for deciding a system of rationality constraints over a 0-cell is then 
an algorithm for deciding whether or not a given real algebraic number P = 
is rational. Naively, one might try to solve this problem in the following way: 

Given P presented as a root-triple {p G Z[a;], / , u), fully factor p over Q[a:]. 

Then, ,9 G Q iff the factorisation of p contains a linear factor of the form 

{x — q) with q G [/, u]. 

From the perspective of theorem proving, the problem with this approach is 
that it is difficult in general to establish the “completeness” of a factorisation. 
While it is easy to verify that the product of a collection of factors equals the 
original polynomial, it can be very challenging (without direct appeal to the 
functional correctness of an implemented factorisation algorithm) to prove that 
a given polynomial is irreducible, i.e., that it cannot be factored any further. 
Indeed, deep results in algebraic number theory are used even to classify the 
irreducible factors of binomials [8]. Moreover, univariate factorisation can be 
computationally expensive, especially when one is only after rational roots. 

We would like the steps in our proofs to be as clear and obvious as possible, 
and to minimise the burden of formalising our procedure as a tactic in a proof 
assistant. Thus, we shall go a different route. To decide whether or not a given a 
is rational, we apply a simple but powerful result from high school mathematics: 

Theorem 4 (Rational roots). Let p{x) = G Z[x] \ {0}. If f G Q 

s.t. p{q) = 0 and gcd{a, b) = 1, then a \ oq and 6 | a„. 

Proof. A straightforward application of Gauss’s lemma. 

Given Theorem 01 we can decide the rationality of a simply by enumerat¬ 
ing potential rational roots qi,... ,qk and checking by evaluation whether any 
qi satisfies {I < qi <r A p{qi) = 0). Then, to decide whether a satisfies a given 
system of rationality constraints, e.g., Pq = [{x^ G Q) A {x ^ Q)], we first com¬ 
pute a root-triple representation for and then test a and for rationality 
as described. This process clearly always terminates. To make this more efficient 
when faced with many potential rational roots, we can combine (i) dividing our 
polynomial p hy {x — q) whenever q is realised to be a rational root, and (ii) 
various cheap irreducibility criteria over Q[a:] for recognising when a polynomial 
has no linear factors over Q[a:] and thus has no rational roots. 


3.2 Deciding integrality constraints 

Integrality constraints over an unbounded 1-cell WLOG let c = ]a,-|-oo[ 
with a > 0. Consider F = /q A Iz with 


' fea 

G Z) A 


. 2=1 


^4 

A(^ 

2=1 


,204(2) 


Z) 


We use the notation (jj : F to mean that the constraint </> is present as a conjunct 
in F. It is convenient to also view F as a set. Let F denote the closure of F 
under the following saturation rules: 


1. (cc" ^Q) ^ (cc" ^Z) -.F 

2. (x" G Z) : C (x" G Q)j C 

3. (x" G Z) : C A (cc'" ^ Z) : C ^ {x ^ <Q) : F_ 

4. (x" G Z) : C A (a:"* G Q) : r ^ (a;™ G Z) : C 

5. (a:" G Z) : C A (a:"* ^ Z) : C (a:™ ^ Q) : F 


This saturation process is clearly finite. The soundness of rules 1 and 2 is obvious. 
The soundness of rules 3-5 is easily verified by the following lemmata. 

Lemma 5 (Soundness: rule 3). (x" G Z) A (x™ ^ Z) —>■ (x ^ Q) 

Proof. Since x™ ^ Z, we know x ^ Z. Suppose x G Q. T/ien x = ^ s.t. 
gcd{a, b) = 1. Thus, a" = x"6". T/ims, 6 | a. Recall gcd{a, 6) = 1. So, 6=1. iJatt 
t6en X = a G Z. Contradiction. □ 


Lemma 6 (Soundness: rule 4). (x" G Z) A (x"* G Q) —>■ (x™ G Z) 

Proof. Let d = deg{x). By Lemma\^ d \ n and d | m. If d = n, then x"* = (x")^ 
for some k £ N and thus x"* G Z. Otherwise, d < n. Let x'^ = f £ Q s.t. 
gcd{a, b) = 1. Thus, x” = (x*^)^ = ^ G Z for some k £ N. But then 6=1, and 

thus x‘^ £ Z. S'o, as d I m, x™ £ h as well. □ 

Lemma 7 (Soundness: rule 5). (x" £ 71) A (x™ ^ Z) —>■ (x™ ^ Q) 

Proof. Assume (x" G Z) and (x™ ^ Z) but (x™ G Q). But then (x™ G Z) 6?/ 

rnZe /(. Contradiction. □ 

Let us now prove that these rule^ are complete for deciding the satisfiabil¬ 
ity of systems of rationality and integrality constraints over unbounded 1-cells. 
Let Fq (resp. Tz) denote the collection of rationality (resp. integrality) con¬ 
straints present in F. Intuitively, we shall exploit the following observation: The 
construction of F projects all information pertaining to the consistency of the 
combined rationality and integrality constraints of F onto Fq. Then, if Tq is 
consistent, i.e., 3d G N satisfying 'D{Fq), this will impose a strict correspondence 
between Fq and Fz. From this correspondence and a least d witnessing T>(Tq), 
we can construct an algebraic real satisfying F. 

^ In fact, the completeness proof shows that rule 3 is logically unnecessary. Neverthe¬ 
less, we find its inclusion in the saturation process useful in practice. 






Lemma 8 (rQ-Fi correspondence). If Fi contains at least one positive in¬ 
tegrality constraint, then 

Vm e N [(x™ e Q) : r ^ (x™ e Z) : T] 

and 

Vm G N [(x™ ^ Q) : r (a;™ ^ Z) : r] . 

Proof. Let us call the first eonjunet A and the seeond B. {A =>) As /z contains 
at least one positive integrality eonstraint, rule 4 guarantees (cc"* € h) : F. {A <=) 
Immediate by rule 2. {B Immediate by rule 1. {B <=) As /z contains at least 
one positive integrality constraint, rule 5 guarantees (cc™ ^ Q) : F. □ 

Theorem 5 (Completeness of T-satnration method). Let F = Fq A F^ 

be a system of rationality and integrality constraints, and c C R an unbounded 
Fcell. Then, T>{Fq) is consistent over N iff F is consistent over c. 

Proof (<^=) Immediate by Theorem\^and the soundness of our saturation rules. 
(=>) We proceed by cases. 

[Case 1: F contains no positive rationality constraint]: Then, by Lemma\^and 
the consistency ofT>[Fq), Iz must contain no positive integrality constraints. 
But then it is consistent with F that every power of x listed in F be irrational. 
Let k G N be the largest power s.t. appears in a constraint in F. Then, by 
Lemma\^ any a G c s.t. deg{a) > k will satisfy F. By Lemma\^ we can always 
find such an a in c, e.g., we can select a G c s.t. deg{a) = fc + 1. 

[Case 2: F eontains a positive rationality constraint but no positive integrality 
constraints]: By the consistency ofTt^Fq), it is consistent with F for every power 
of X listed in F to be non-integral. Let d G^ be the least natural number satisfying 
'D{Fiq). Then, we can satisfy F with an a s.t. deg{a) = d with G (Q \ Z) for 
each x’^^ appearing in a constraint in F. By Lemma\^ we know such an a is 
present in c of the form a = for primes p g. 

[Case 3: F contains both positive rationality and integrality constraints] By 
Lemma\^ the rows of Fq and Tz are in perfect correspondence. Let d G N 6e the 
least natural number satisfying T>{Fq). Since Fq is consistent, we can satisfy F 
by finding an a G c s.t. G Z for every x‘^^ appearing in a constraint in F. 
Recall c is unbounded towards +oo. Thus, c contains infinitely many primes p 
s.t. f/p G c. Let p G c be such a prime. Then, x‘^ — p G Q[x] is irreducible by 
Eisenstein’s criterion. Thus, ffp G c and satisfies F. □ 

Integrality constraints over a bonnded 1-cell Let us now consider the 
satisfiability of T = Jq A Jz over a bounded 1-cell c C K. Given the results of 
the last section, it is easy to see that if 'D{Fq) is unsatisfiable over N, then F is 
unsatisfiable over c. However, as F is bounded on both sides, it is possible for 
'D{Fq) to be satisfiable over N while F is unsatisfiable over c. That is, provided 
T>{Fq) is consistent over N, we must find a way to determine if c actually contains 
some a s.t. F{a) holds. Afterall, even with Fq satisfied over c, it is possible that 
c itself is not “wide enough” to satisfy the integrality constraints Tz. 


WLOG, let c = ]a,/3[ s.t. 0 < a < P € Ko/g- Let 25 (Gq) be satisfied by d G N. 
If r contains no positive integrality constraints, then we can reason as we did in 
the proof of Theorem [5] to show F is satisfied over c. The difficulty arises when 
a positive constraint G Z) appears in Iz. We can solve this case as follows. 

Theorem 6 (Satisfiability over a bounded 1-cell). Let Iz contain at least 
one positive integrality constraint. Let T>(Tq) be satisfiable over N with d G N 
the least witness. Let c = ]q;,/3[ s.t. 0 < a < /3 G Maig. Then, F is satisfiable 
over c iff 3 zG fl Z) s.t. x‘^ — z £ Z[x] is irreducible over Q[x]. 

Proof. (=>) Assume F is satisfied by a £ c. Then, by soundness of F saturation, 
F is satisfied by a as well. By Lemma 0 G Z) : T. Moreover, d is the least 
natural number with this property. As Q < a < P, {r^ | r G c} = \a^,P‘^[. Thus, 
as F is satisfied by a £ c, there must exist an integer z £ ]a‘^,P'^[ s.t. deg{y/z) = 
d. But then by unigueness of minimal polynomials, x‘^ — z is irreducible over 
Q[x]. (<^=) Assume z £ fl Z) s.t. x‘^ — z is irreducible over Q[x]. Let 

7 = ffz and note that 7 G ]a, P[. By Lemmo0 deg(j) = d. Thus, Fq is satisfied 
by 7 . As 7 ^^ G Z, it follows by Lemmo0 that F is satisfied by 7 as well. □ 

By Eisenstein’s criterion, we obtain a useful corollary. 

Corollary 2. Let L>(Tq) be satisfiable with d £ N the least natural number 
witness. Let c = ]a,P[ s.t. Q < a < P £ M.aig- Then, F is satisfiable over c if 
3p £ ]a‘^,P‘^[ s.t. p is prime. 

These results give us a simple algorithm to decide satisfiability of F over c: 
If 21 (Tq) is unsatisfiable over N, then F is unsatisfiable. Otherwise, let d G N be 
the minimal solution to 'D{Fq). Gather all integers {zi, ..., Zk} in J = P'^[. 

If any Zi is prime, F is satisfied over c. Otherwise, for each Zi, form the real 
algebraic number ffzl and check by evaluation if it satisfies F. By Theorem [ 6 l 
F is satisfiable over c iff one of the ffzi £ c satisfies this process. 

Integrality constraints over a 0-cell Finally, we consider the case of T = 
Fq a Fz over a 0-cell [a]. Clearly, F is satisfied over c iff T is satisfied at a. 
By the soundness of T-saturation, if 'D(Fq) is unsatisfiable over N, then F is 
unsatisfiable over c. Thus, we first form F and check satisfiability of 'D{Fq) over 
N. Provided it is satisfiable, we then check F(x ^ a) by evaluation. 

4 Examples 

We have implementec^l our decision method in a special version of the Meti- 
Tarski theorem prover m- We do not use any of the proof search mechanisms 
of MetiTarski, but rather its parsing and first-order formula data structures. 

In the examples that follow, all output (including the prose and format¬ 

ting) has been generated automatically by our implementation of the method. 

^ The implementation of our procedure, including computations over r-cells, F- 
saturation and the proof output routines can be found in the RCF/ modules in the 
MetiTarski source code at http://metitarski.googlecode.com/ 



4.1 Example 1 

Let us decide 3x{(p{x) A r{x)), where 

(/? = — 2 = 0) and F = (x G Q). 

We first compute F, the closure of F under the saturation rules: 

r = (x G Q). 

Observe 21 (Eq) is satisfied (minimally) by c? = 1. 

We next compute an r-cell decomposition of M induced by ip, yielding: 

1. ]-oo, i?oot(x2 - 2, [-2,-l/3])[, 

2. [i?oot(x2 - 2, [-2,-l/3])], 

3. ]Root{x^ - 2, [-2,-l/3]),0[, 

4. [0], 

5. ]0,i?oot(x2 - 2, [l/3,2])[, 

6 . [i?oot(x2 - 2, [1/3,2])], 

7. ]i?oot(x^ — 2, [1/3, 2]),+oo[. 

By IVT, (fi has constant truth value over each such r-cell. Only two r-cells in the 
decomposition satisfy (p: 

[Rootix^ - 2, [-2, -1/3])], [i?oot(x2 - 2, [1/3, 2])]. 

Let us now see if any of these r-cells satisfy F. 

1. We check if [i?oot(x^ — 2, [—2, —1/3])] satisfies F. 

(a) Evaluating (a G Q) for a = Root{x‘^ — 2, [—2, —1/3]). We shall determine 
the numerical type of a. Let p(x) = x^ — 2. By RRT and the root interval, 
we reduce the set of possible rational values for a to {—1, —2}. But none 
of these are roots of p{x). Thus, a G (R \ Q). 

So, the r-cell does not satisfy F. 

2. We check if [Root{x‘^ — 2, [1/3, 2])] satisfies F. 

(a) Evaluating (a G Q) for a = Root{x'^ — 2, [1/3, 2]). We shall determine the 
numerical type of a. Let p(x) = x^ — 2. By RRT and the root interval, 
we reduce the set of possible rational values for a to {1,2}. But none of 
these are roots of p{x). Thus, a G (R \ Q). 

So, the r-cell does not satisfy F. 

Thus, as all r-cells have been ruled out, the conjecture is false. □ 

4.2 Example 2 

Let us decide 3x{p{x) A F{x)), where 

p = True and F = (x^ &%) f\ (x® ^ Z) A (x G Q). 

We first compute T, the closure of F under the saturation rules: 

F = {x ^ Z)A(x G Z)A(x ^ Q)A(x G Q)A(x^ G Q)A(x^ G Z)A(x^ ^ Q)A(x® ^ Z). 


But, F is obviously inconsistent. Thus, the conjecture is false. 


□ 



4.3 Example 3 

Let us decide 3x{(p{x) A r{x)), where 

ip = ((x^ — 7 > 3) A + X + 1 < 50)) and F = (x^ ^ Q) A (x^ € Z). 

We first compute F, the closure of F under the saturation rules: 

F = (x^ ^Q)A (x^ ^ Z) A (x^ G Z) A (x^ G Q). 

Observe 'D{Fq) is satisfied (minimally) by c? = 3. 

We next compute an r-cell decomposition of M induced by p, yielding: 

1. ]—oo, Root{x^ + X — 49, [—8, —1/50])[, 

2. [i?oot(x2+x-49, [-8,-l/50])], 

3. ]i?oot(x2+x-49, [-8,-l/50]),0[, 

4. [0], 

5. ]0,i?oot(x3-10, [57/44,5/2])[, 

6 . [Root{x^ - 10, [57/44,5/2])], 

7. ]Root{x^ - 10, [57/44,5/2]),i?oot(x2 + x-49, [401/100,8])[, 

8 . [i?oot(x2+x-49, [401/100,8])], 

9. ]Root{x^ + X - 49, [401/100,8]), +oo[. 

By IVT, p has constant truth value over each such r-cell. Only one r-cell in the 
decomposition satisfies p: 

]Root{x^ - 10, [57/44, 5/2]), Root{x"^ -f x - 49, [401/100, 8])[. 

Let us now see if any of these r-cells satisfy F. 

1. We check if ]Root{x^ - 10, [57/44,5/2]), Rootix"^ -f x - 49, [401/100,8])[ sat¬ 
isfies F. Call the boundaries of this r-cell L and U. As F contains a positive 
integrality constraint and d = 3, any satisfying witness in this r-cell must be 
of the form for z an integer in ]L^,U^[. The set of integers in question 
is Z = {z G Z I 11 < z < 276}, containing 266 members. We shall examine 
^ for each z G Z in turn. 

(a) Evaluating (a^ ^ Q) for a = Root{x^ — 11, [1/12,11]). Observe = 
Root{x^ — 121, [1/144,121]). We shall determine the numerical type of 
a^. Let p(x) = x^ — 121. By RRT and the root interval, we reduce the 
set of possible rational values for to {1,11,121}. But none of these 
are roots of p{x). Thus, o? G (R \ Q). 

(b) Evaluating [c? G Z) for a = Root{x^ — 11, [1/12,11]). Observe = 
Root{x^ — 1331, [1/1728,1331]). We shall determine the numerical type 
of . Let p(x) = x^ — 1331. By RRT and the root interval, we reduce 
the set of possible rational values for to {1,11,121,1331}. Thus, we 
see = 11 G Z. 

Witness found: Root{x^ — 11, [1/12,11]). So, the r-cell does satisfy F. 


Thus, the conjecture is true. 


□ 



5 Discussion and Related Work 


Let us describe some related results that help put our work into context. 

— The existence of rational or integer solutions to univariate polynomial equa¬ 
tions over Q[a;] has long been known to be decidable. The best known algo¬ 
rithms are based on univariate factorisation via lattice reduction . 

— Due to Weispfenning, the theory of linear, multivariate mixed real-integer 
arithmetic is known to be decidable and admit quantifier elimination [20) . 

— Due to van den Dries, the theory of real closed fields extended with a 
predicate for powers of two is known to be decidable Avigad and Yin 
have given a syntactic decidability proof for this theory, establishing a non¬ 
elementary upper bound for eliminating a block of quantifiers [5]. 

— Due to Davis, Putnam, Robinson and Matiyasevich, the 3^ nonlinear, equa- 
tional theories of arithmetic over N and Z are known to be undecidable 
(“Hilbert’s Tenth Problem” and reductions of its negative solution) [11) . 

— The decidability of the 3^ nonlinear, equational theories of arithmetic over 
N and Z is open. 

— Due to Poonen, the V^3^ theory of nonlinear arithmetic over Q is known 
to be undecidable m- This is an improvement of Julia Robinson’s original 
undecidability proof of Th{Q) via a V^3’^V® definition of Z over Q [IB] . 

— Due to Koenigsmann, the and theories of nonlinear arithmetic 

over Q are known to be undecidable, via explicit definitions of Z over Q [mo]. 

— The decidability of the 3* equational nonlinear theory of arithmetic over Q 
is open for fc > 1 (“Hilbert’s Tenth Problem over Q”). 

Our present result — the decidability of the nonlinear, univariate theory of the 
reals extended with predicates for rational and integer powers — fills a gap 
somewhere between the positive result on linear, multivariate mixed real-integer 
arithmetic, and the negative result for Hilbert’s Tenth Problem in three variables. 

Next, we would like to turn our decision method into a verified proof proce¬ 
dure within a proof assistant. The deepest result needed is the Prime Number 
Theorem (PNT). As Avigad et al have formalised a proof of PNT within Is- 
abelle/HOL [T], we are hopeful that a verified version of our procedure can be 
built in Isabelle/HOL [14] in the near future. To this end, it is useful to observe 
that PNT is not needed by the restriction of our method to deciding the ratio¬ 
nality of real algebraic numbers like and -x/S -I- \/5. Thus, a simpler tactic 
could be constructed for this fragment. 

Finally, we hope to extend the method to allow constraints of the form 
{p{x) G Q) for more general polynomials p{x) G Z[x]. The key difficulty lies 
with Lemma [2] This crucial property relating the degree of an algebraic number 
to the rationality of its powers applies to “binomial root” algebraic numbers, 
but not to algebraic numbers in general. For example, consider a = ^2 + \f2. 
Then, the minimal polynomial of a over Q[x] is x^ — 4x^ — 8x + 2, but ^ Q. 
Thus, in the presence of richer forms of rationality and integrality constraints, 
our degree constraint reasoning is no longer sufficient. We expect to need more 
powerful tools from algebraic number theory to extend the method in this way. 




6 Conclusion 


We have established decidability of univariate real algebra extended with pred¬ 
icates for rational and integer powers. Our decision procedure combines com¬ 
putations over real algebraic cells with the rational root theorem and results 
on the density of real algebraic numbers. We have implemented the method, 
instrumenting it to produce readable proofs. In the future, we hope to extend 
our result to richer systems of rationality and integrality constraints, and to 
construct a verified version of the procedure within a proof assistant. 

Acknowledgements. We thank Jeremy Avigad, Wenda Li, Larry Paulson, Andras 
Salamon and the anonymous referees for their helpful comments. 
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